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Recent high-redshift Type la supernovae results can be used to set new bounds on a possible vari- 
ation of the gravitational constant G. If the local value of G at the space-time location of distant 
supernovae is different, it would change both the kinetic energy release and the amount of 56 Ni 
synthesized in the supernova outburst. Both effects are related to a change in the Chandrasekhar 
mass Mc\i <x G~ 3//2 . In addition, the integrated variation of G with time would also affect the 
cosmic evolution and therefore the luminosity distance relation. We show that the later effect in the 
magnitudes of Type la supernovae is typically several times smaller than the change produced by 
the corresponding variation of the Chandrasekhar mass. We investigate in a consistent way how a 
varying G could modify the Hubble diagram of Type la supernovae and how these results can be used 
to set upper bounds to a hypothetical variation of G. We find G/Go < 1.1 and G/G < 10~ 11 yr~ 1 
at redshifts z ~ 0.5. These new bounds extend the currently available constrains on the evolution 
of G all the way from solar and stellar distances to typical scales of Gpc/Gyr, i.e. by more than 15 
orders of magnitudes in time and distance. 



PACS numbers: 98.80.Cq, 04.50.+h 
I. INTRODUCTION 

One of the most important challenges of modern 
physics is the quantization of the gravitational force. 
The undergoing attempts to create such theories has re- 
opened the subject of varying fundamental constants. To 
this regard it is worth noticing that the constancy of the 
fundamental constants, and of the gravitational constant 
in particular, has been questioned for a long time 
and that early attempts to unify gravity with electromag- 
netism predicted such kind of variations. Although 
modern theories, like the string theory and the M-thcory 
(see |y for a recent review), do not necessarily require 
a variation of the fundamental constants they provide a 
natural and self-consistent framework for such variations 
(see Q and jl0| for excellent descriptions of the theo- 
retical background). As a general result, modern theo- 
ries predict that in the ordinary three-dimensional sub- 
space, gauge couplings like the fine structure constant a 
or the gravitational constant G should vary as the inverse 
square of the mean scale of the extra dimensions. Hence, 
the evolution of the scale size of the additional dimen- 
sions is related to the variation of fundamental constants 
|TT| |i"3[ . Moreover it has been recently shown jl3| that 



a cosmological variation of a may proceed at different 
rates at different locations in space-time. The way in 
which the time variations of a and G are linked is model 
dependent but a typical relation is: Aa/a 2 ~ AG/G. 
There have been several attempts to measure the rate 
of variation of a, providing different results for different 
look-back times. For instance, |jl4j used the recently re- 
leased CMB anisotropy data to set up early-universe con- 
straints on a time-varying a and found no evidence for 
such a change, whereas |n]|n| used high resolution spec- 
troscopy of QSO absorption systems to find statistical ev- 
idences for a smaller a at a redshift range 0.5 < z < 3.5. 
Of course, since a cosmological variation of a (and, conse- 
quently, of G) can proceed at different rates for different 
redshifts jl3| both studies are not necessarily in conflict. 
There have been also many attempts to measure a time 
variation of the gravitational constant which will be dis- 
cussed later in §V. For the moment it is important to 
mention here that most of these bounds come either from 
local measurements (the sun, our solar system or the so- 
lar neighborhood) or from very early times measurements 
(namely Big-Bang nucleosynthesis), whereas at interme- 
diate look-back times there are not such measurements. 
Type la supernovae (SNIa) are supposed to be one of 
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the best examples of standard candles. This is because, 
although the nature of their progenitors and the detailed 
mechanism of explosion are still the subject of a strong 
debate, their observational light curves are relatively well 
understood and their individual intrinsic differences can 
be accounted for. Under these assumptions, thermonu- 
clear supernovae are well suited objects to study the Uni- 
verse at large, especially at high redshifts (z ~ 1.0), 
where the rest of standard candles fail in deriving reli- 
able distances, thus providing an unique tool for deter- 
mining cosmological parameters or discriminating among 
different alternative cosmological theories. 

Using the observations of high redshift Type la su- 
pernovae (z > 0.1) and low redshift (z < 0.1) super- 
novae, both the Supernova Cosmology Project (l(| and 
the High-z Supernova Search Team jl7| found that the 
peak luminosities of distant supernovae appear to be 
~ 0.20 magnitude fainter than predicted for an empty 
universe and 0.25 magnitude fainter than predicted by 
a standard decelerating universe, with a presumed mass 
density Qm — 0.3. To be more precise, at the la con- 
fidence level, the results of both groups can be well ap- 
proximated by the relation: 



o.8 n M ~0.6fl A = -0.2 ±0.1. 



(1) 



However these conclusions rely on the assumption that 
there is no mechanism able to produce an evolution of 
the observed light curves over cosmological distances. In 
other words: both teams assumed that the relation be- 
tween the intrinsic peak luminosity and the time scales of 
the light curve were exactly the same for both the low-z 
and the high-z supernovae. The possible consequences 
for evolutionary effects in SNIa due to changes in the 
zero age mass and metalicity of the progenitor star have 
been explored by several authors ||l8|-^o| , who found that 
changes in the underlying population cause a change in 
the maximum brightness by about 0.1-0.2 magnitudes. 

The SNIa results have already motivated a significant 
number of papers that search for bounds on the vari- 
ation in fundamental constants |2l|-|25|] or new cosmo- 
logical scenarios, such as quintessence models [^6| and 
scalar- field cosmologies [^7j . This burst of interest is due 
to the conceptual problems that arise from infering the 
existence of a cosmological constant A ~ 10~ 122 c 3 /G7?t 
or facing the cosmic (dark) matter problem (see pq] and 
references therein). There have been many suggestions 
that the apparent complications that arise can be elimi- 
nated by modifying the laws of gravity |2j]-|3(| . 

Recent cosmological observations, such as the lastest 
CMB Boomerang and Maxima data |B7 [ 
flat Universe: flu — 0, i.e. Qm + = 
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indicate a 
This result, 



together with the above equation ^ points in the direction 
of a non-cero A, although other interpretations are also 
possible HQ. 

The purpose of this paper is to analyze the effect of 
varying G in the current interpretation of the Hubble di- 



agram of distant SNIa and to use this analysis to set up- 
per bounds on its rate of change. The paper is organized 
as follows: in §11 we describe the effects of a varying G on 
the physics of supernovae; in §111 and in the Appendix, 
we analyze the effects of a varying G on the luminosity 
distance of distant supernovae. In §IV we present a like- 
lihood analysis of the SNIa data which is then used to set 
upper bounds in the evolution of G. Finally in section V 
we discuss our results and draw our conclusions. 



II. THE EFFECTS OF A VARYING G ON THE 
PHYSICS OF SUPERNOVAE 

Simple analytical models of light curve (see, for in- 
stance, ]4l|) predict that the peak luminosity is propor- 
tional to the mass of nickel synthesized, which in turn, 
to a good approximation, is a fixed fraction of the Chan- 
drasekhar mass {Mm oc Mch), which depends on the 
value of gravitational constant as: Mch oc G~ 3 / 2 . The 
actual fraction varies when different specific SNIa scenar- 
ios are considered (e.g. [^2|,^3)), but the physical mech- 
anisms relevant for type la supernovae naturally relates 
the energy yield to the Chandrasekhar mass. Here we 
will only focus on Chandrasekhar mass models since there 
are growing evidences that the sub-Chandrasekhar mass 
models do not fit well the observations (see, for instance, 
p4| ). In summary, whatever the actual scenario is, we 
will assume that the same mechanism for the ignition and 
the propagation of the burning front is valid for SNIa at 
high and low redshifts. Thus, since the peak luminosity is 
proportional to the total amount of nickel synthesized in 
the supernova outburst we have L oc G~ 3 / 2 , and, there- 
fore, for a slow decrease of G with time, distant super- 
novae should be dimmer than predicted for a standard 
scenario. Under this assumptions, we have: 



15, { G \ 
M-M = T log(-) 



(2) 



where M stands, as usual, for the absolute magnitude, 
G for the precise value of the gravitational constant at 
a given redshift and the subscript denotes their local 
values. Note that the above equation does not require 
knowledge of the (unknown) proportionality constant re- 
lating the supernova luminosity with the Chandrasekhar 
mass. This dependence factorizes out under the assump- 
tions above, and the final differential result is only sen- 
sitive to the values of G. From this equation we can 
see that in order to reduce the apparent luminosity of 
distant supernovae by Am ~ 0.2 a dramatic change of 
G is required: G/Gq ~ 1.13. This value should be re- 
garded as an upper bound in the sense that part or all 
of the Am ~ 0.2 difference found by Jl6| . |T7| | could be at- 
tributable to the hypothesis of an accelerating universe. 

In order to test the validity of our argument we have 
computed a series of models of Type la supernovae ex- 
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plosion and their corresponding light curves, according 
to the procedure described in J45[ and references therein, 
with the present local value of G, 1.1 and 1.2 times this 
value. The explosion model was a delayed detonation 
starting from a central density of 2.0 x 10 9 g/cm 3 a core 
temperature of 2.0 x 10 8 K and making the transition 
from deflagration to detonation when the flame density 
went below 2.0 x 10 7 g/cm 3 . Our study is based on de- 
layed detonation models, because these have been found 
to reproduce the optical and infrared light curves and 
spectra of Type la supernovae reasonably well |l^,^ P| ■ 
The model parameters, ignition density and transition 
density, are those that allow us to reproduce a typical 
Type la Supernovae. The results are shown in Table 1 
and Figure 1. 
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FIG. 1. Bolometric light curves of SNIa for the local value 
of the gravitational constant, Go (solid line), for Gi = l.lGo 
(dotted line) and for Gi shifted upwards by 0.18 magnitudes. 
See text for details. 

In Table 1 we show the mass of the white dwarf model 
in hydrostatic equilibrium from which the explosion was 
computed, M-yvDj the kinetic energy, K, the mass of syn- 
thesized nickel, M-^i, the peak bolometric magnitude, 
Muff , and a measure of the width of the light curve, 
Amis — see below for a precise definition. All the masses 
are expressed in solar units. In Figure 1 we show the light 
curves for G/G = 1.0 — solid line — G x = G/G = 1.1 
— dotted line — and this last light curve shifted upwards 
by a constant amount of 0.18 magnitudes — dashed- 
dotted line — in accordance with the behavior predicted 
by Eq.|]. 

As it can be seen in Table 1 the simple energetic ar- 
gument presented above is in good agreement with the 



detailed calculations presented here and, thus, the en- 
ergy liberated in the supernova outburst indeed scales as 
G~ 3 / 2 . Hence, the bolometric magnitude at the peak of 
the light curve computed with a larger value of G should 
be moved upwards by a fixed constant amount which de- 
pends on the exact value of G. Figure 1 clearly shows 
this behavior. It should be stressed at this point that 
our analysis does not depend on the detailed physics of 
Type la supernovae since the functional dependence on 
the Chandrasekhar mass comes from basic physical ar- 
guments and it is unavoidable, no matter which are the 
(unknown) details of the explosion unless we change the 
physics underlaying the Chandrasekhar mechanism. 

Figure 1 also shows that once this vertical shift is done, 
the duration of the supernova outburst is also modified 
by a varying G, being the decline faster for the models 
with larger G, especially at late times. As it can be seen 
there, in the region near the maximum (see the insert for 
a close up of this region) the difference between both light 
curves is small when due account of the vertical shift is 
done. A good parameterization of the slope (and, thus, 
of the width) of the light curve is Amis, which is defined 
as the difference in the apparent magnitude 15 days af- 
ter the maximum. It turns out that the maximum of 
the light curve for the models presented in Fig. 1 occurs 
at 13.0 and 12.5 days, being Ami 5 = 0.85 and 0.89, re- 
spectively. This, in turn, implies that since the template 
light curve used to calibrate the distances to distant su- 
pernovae takes into account the width of the light curve, 
and in particular Amis is used, a variation of the gravi- 
tational constant should ultimately affect as well the final 
value of the derived distances. However, given the slight 
variation of the Amis parameter this can be clearly con- 
sidered as a second order effect. Nevertheless, as it can 
be clearly seen in Fig. 1, although there is not any large 
difference in the rise-times or in the early time light curve 
due to a varying G there is indeed an appreciable differ- 
ence in the overall duration of the supernova event. To 
be precise the widths at Mboi = — 16 are 126 and 112 
days, respectively. 

This result can be interpreted in terms of very simple 
physical considerations, and in particular, in terms of the 
simplified model of light curve presented in fli"|| which has 
a reasonable accuracy (of the order of ~ 20%). According 
to this analytic model of light curve, the width of the 
peak of the light curve of SNIa is given by: 

TABLE I. Overall characteristics of the supernova explo- 
sion for several values of G, and the same ignition density. 
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Mt, \ 1/4 



(3) 



where M c j is the ejected mass and K is the kinetic en- 
ergy. Within our current knowledge of the mechanisms 
of explosion of SNIa both the ejected mass and the super- 
nova kinetic energy can be considered proportional to the 

1/2 

Chandrasekhar mass, and therefore we have r oc M Ch 
or, equivalently, r oc G~ 3 / 4 . Thus we have 
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The ratio of the durations of the supernova outburst 
calculations presented above matches reasonably well the 
behavior predicted by Eq. Hence, in the case in which 
a varying gravitational constant is considered, the overall 
time scales of the supernovae light curves should depend 
as well on the actual value of G. 

It has been recently claimed that there is a mean evo- 
lution in the rise-times of local and distant supernovae 
|E3,p3l. In particular the widths of the light curve when 
the supernova is 2.5 magnitudes fainter than the peak 
luminosity was found to be ro = 45.0 ± 0.15 (at z ~ 0) 
and r = 43.8 ± 0.40 days (at z ~ 0.5), were the er- 
rors in the widths were ascribed solely to the errors in 
the rise-times. Using this data and Eq.Q, pj| obtained 
G/G < 1.037 ± 0.017 (2cr errors) at z - 0.5, a vari- 
ation of G very similar to the one needed to explain 
the change in the peak luminosity. Subsequent studies 
p4j have demonstrated that the rise-times of local and 
distant supernovae are consistent each another and that 
the rise-time uncertainties were underestimated, being 
revised upwards to ±1.2 days statistical and ti.g days 
due to systematical bias under extreme situations. Ac- 
cording to this last analysis there is still some room for 
deviations of the light curve of high redshift supernovae 
at late times. This is exactly what we have found. How- 
ever, it is worth mentioning at this point that according 
to this study these late time deviations systematically in- 
fluence the rise-time determinations. Conversely, should 
we have tried to fit the late time light curve we would 
have found a significant difference in the corresponding 
rise-times. In the light of this new analysis, there is no 
significant evidence for a possible change in G. For our 
purposes, since the situation is not yet clear from the ob- 
servational point of view our analysis will only rely on 
the limits imposed by Eq.Q. Note however that, on the 
other hand, should a firm estimate of the maximum value 
of the difference between the local and distant supernovae 
timescales could be eventually obtained a very stringent 
upper limit to the rate of variation of the gravitational 
constant would be derived, which, additionally, would 
not depend on the adopted cosmological model but on 
the well calibrated relation between the duration of the 
supernova outburst and its peak magnitude. 



III. THE HUBBLE DIAGRAM 

The next question we want to address here is how an 
hypothetical variation of the gravitational constant (of 
the amount of a few percent) translates in the Hubble 
diagram of distant supernovae. In addition to the change 
in the intrinsic energy release and of the duration of the 
supernova event induced by the variation in the Chan- 
drasekhar mass, the integrated evolution of G with time 
would also affect the cosmic evolution and, therefore, the 
luminosity distance relation. To make this quantitative 
we need to consider non-standard cosmological models. 
It is however beyond the scope of this paper to discuss 
specific theoretical models to replace the standard theory 
of General Relativity. This has been discussed in detail 
elsewhere in the literature ||55|-p0|, and also in the context 
of an accelerating universe [p5l6lH6J 




FIG. 2. Hubble diagrams for the high-redshift SNIa data 
with different choices of cosmological parameters and different 
values of G . Continuous lines correspond to the standard 
cosmology with G(z) — Go. Dashed lines bracket the effects 
on the peak luminosity from a ~ 5% variation in G. The 
observational data has been obtained from [ p^[ . 

In the appendix we show in a self-consistent way that 
for plausible models which incorporate a varying G, such 
as scalar-tensor theories (STTs), the possible effect of a 
varying G on the cosmological evolution gives a contribu- 
tion to the luminosity distance relation at the distances of 
interest which is several times smaller than the effect pro- 
duced by the same variation of G on the Chandrasekhar 
mass. We will therefore concentrate our analysis on the 
effects of the physics of the supernovae. This comple- 
ments the analysis of [[35 62 6^] which neglected the ef- 
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fects in the physics of the supernovae and used the evi- 
dence for the accelerating universe as a way to constraint 
cosmic evolution in non-standard theories of gravity. 

In analogy to STTs, we will parameterize the evolution 
of G in terms of the strength of coupling parameter, uj, 
as: 



G(z) = G (l + 2)^+^17, 



(5) 



which provides our definition for uj(z) — see the appendix 
for further details. Thus a value of uj ~ 10 produces a 
~ 4% increase in G at z ~ 0.5, while u> ~ — 10 produces a 
~ 5% decrease in G at z ~ 0.5. In summary, Eq.J| gives 
the change in G as a function of uj while Eq.[A3-A5] give 
the corresponding cosmic evolution. Thus we have: 



m(z) = M + 5 log d L + 25 
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4(1 + u) 



log(l 



(6) 



where dt — d^(z, Hm, ^a) is obtained from the line-of- 
sight comoving distance (see the appendix) . For illustra- 
tive purposes Figure || shows the above relation for two 
representative cosmological models, including the effects 
of uj in rfi, for uj = ±10 (dotted lines), which correspond 
to a change of G of ~ 5%, and the standard (uj = oo) 
case (solid line). 



point that we have also tried linear fits to to(z) ~ ujq+uj' z 
and found equivalent results. 
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FIG. 3. Confidence contours for the best fit SNIa data in 
the plane (G/Go, ^a) for a flat universe f2.R = 0. 



IV. BOUNDS ON G/G 

We have re-done the likelihood analysis of the Su- 
pernova Cosmology Project allowing for a varying G. 
Thus, we use the same observational data but with the 
magnitude-distance relation given by Eq.^j. Here we 
have one extra function uj(z) to be fitted. At low red- 
shifts the last term in the r.h.s. of Eq.[|| has a negligible 
contribution. Given that most of the SNIa at high red- 
shift cluster around z ~ 0.5 and that we are expecting 
uj to be large (so that G ~ Go in Eq.||), the effect of 
uj(z) in the fit to Eq.|| is dominated by the value of uj 
at the mean redshift of the SNIa sample, z ~ 0.5. Hence, 
we can approximate uj(z) ~ uj(z ~ 0.5) and fit the data 
as a function of this new extra parameter, uj(z ~ 0.5) or 
G(z~0.5)/G(z~0) = G/G . 

Our results are plotted in Figures 3 and 4, where we 
show the confidence contours (at the 99%, 90%, 68% 
- solid lines — 5% and 1% confidence level — dot- 
ted lines) obtained from the fit to the Hubble diagram 
of SNIa. Figure 3 shows the likelihood contours for 
G/G = G(z ~ 0.5)/G(z ~ 0) as a function of fl A 
for a flat £Ir — universe, whereas Figure 4 shows 
the likelihood contours in the (G /Gq,VLm) plane for the 
case = 0. As it can seen in these figures, the ex- 
pected departures from the standard case (G/Gq ~ 1) 
are quite small for a reasonable choice of cosmological 
parameters. This, in turn, justifies our approximation 
w(z)~w(z~0.5). It is also worth mentioning at this 



n A =o 




FIG. 4. Confidence contours for the best fit SNIa data in 
the plane (G/Gq,Q.m) for £2a = 0. 

The confidence intervals of these figures can be used 
to set bounds in G/Go- These are bounds in the sense 
that, given a cosmological model, we assume that all the 
difference in SNIa corrected peak luminosities can be at- 
tributed solely to a difference in Chandrasekhar mass. 
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For example, for the flat Q A 
confidence level: 



model we have at la 



Gq ~ 



n A ~ o.o , n 



M 



while for the flat il A — 0.8 case 
— < 1.08 ; n A ~ 0.8 



Af 



1.0, 



0.2 



(7) 



(8) 



also at la confidence level. In terms of the u> ~ u(z ~ 
0.5) parameter these later bounds translate into cj > 1.2 
for 17a — and cj > 4.3 for Q\ ~ 0.8. It is important 
to mention here that these constraints are quite loose 
when compared with the bounds on oj from very long 
baseline interferometry in the local Universe, ujq > 2600, 
fl60| , but, on the other hand, are new in the sense that 
they correspond to an earlier (or equivalently more dis- 
tant) Universe. Moreover, the SNIa data provide more 
interesting limits to G/G. To obtain them we can use 



H'-f)<">- 

where At is the look-back time to z ~ 0.5: 

dz' 

At = 



(9) 



(10) 



with H(z) given by Eq. |A3|| . For a flat Q A ~ model 
we have that At ~ 4.6 x 10 9 yr//i7o (hro is the Hubble 
constant in units of Hq = 70 Km/s/Mpc), while for the 
flat f^A ~ 0.8 case we have At ~ 6.0 x 10 9 yr//i7 . Thus 
we find: 



(l + z>)H(z>) 



< 36 x 10" 12 h 70 /yv ; ft A ~ 0.0 , fi M ~ 1.0 (11) 



G 
G 

- < 12 x 10~ 12 h 70 /yr ; fl A 0.8 , Q M 0.2 (12) 

It should be stressed that these are bounds on G/G 
around z ~ 0.5. Several local bounds on the rate of 
change of the gravitational constant, Go/Go, have been 
obtained, for example, from binary pulsars, from the 
Viking Radar and from Lunar Laser Ranging, yielding 
typical upper bounds of Go /Go < 10 x 10~ 12 yr _1 (see, 
for instance, ]60[| ). Other bounds come from the theory 
of stellar evolution, like white dwarf cooling |p4{| , being 
the bounds of the order of 10 x 10~ 12 yr _1 . Finally, it 
should be mentioned that some other local bounds are as 
low as Go /Go < 6 x 10~ 12 yr _1 |65|. Note, however, that 
the values of all these bounds are comparable to that ob- 
tained here. Thus, despite the fact that the SNIa bounds 
on G/Go are quite loose, the longer time baseline ob- 
tained by using high redshift measurements puts stronger 
constraints on G/G. These bounds on the change in G 
correspond to redshifts which have not been explored yet 
in the more standard gravitational tests, mostly based 



on solar system and stellar physics. We can combine the 
SNIa and local bounds further to set a bound on G: 



| < 4 x 10~ 21 (Wyr) 2 



(13) 



In terms of lo(z) we can also set some further limits by 
doing a Taylor expansion as in J62[ : 



+ z 



cLj 



dz 



so that we find: 



duj~ 



dz 
duj- 1 



< i.7 ; n A ~ o.o , n 



M 



dz 



< 0.5 ; n A ~ 0.8 , n 



M 



1.0 



0.2 



(14) 

(15) 
(16) 



Future experiments, such as the Supernovae Accelera- 
tion Probe (SNAP), could achieve a few percent magni- 
tude errors up to redshifts of z ~ 1.7 (see f56|). By then, 
we can also fairly assume that other observational data, 
such as the LSS and CMB experiments, will provide a 
knowledge of a few percent on the cosmological param- 
eters [Eg] , Thus we can translate the uncertainty in the 
magnitude directly into bounds on G: 



Am 



15 
T 



(17) 



Note that the possible effects of a varying G on the 
luminosity distance discussed in the Appendix, are still 
several times smaller than the above contribution for z ~ 
1.7. For example, a 1% uncertainty in Am (both from 
peak luminosity and cosmological parameter errors) will 
give us a 0.6% bound in G/Go < 1.006 at z — 1.7 or, 
equivalently, uj > 167. The later value is not particularly 
impressive as compared to the local bounds in the context 
of Jordan-Brans-Dicke (JBD) theories (u> > 2500). But 
note that z ~ 1.7 corresponds to a look-back time in 
Eq.0 as large as At ~ 10 x 10 9 yr//i 70 . Thus, future 
data would eventually yield firm bounds for G as low as: 



^ < 6 x 10- 



/Wyr 



1.7 



(18) 



which are more competitive than the current local values. 
The bounds to the change in cj, could be reduced to: 



duj 



dz 



< 3.5 x 10" 



1.7 



(19) 



Finally it is interesting to mention here that our anal- 
ysis has been restricted to the peak magnitudes of su- 
pernovae. Comparable bounds can be found from the 
duration of supernovae events (Eq.Q), if the statistical 
and systematic errors are reduced significantly, being the 
advantage of these last ones that are independent of the 
adopted cosmological model. 
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V. DISCUSSION AND CONCLUSIONS 

In astrophysics and cosmology the laws of physics (and 
in particular the simplest version of general relativity) are 
extrapolated outside its observational range of validity. It 
is therefore important to test for deviations of these laws 
at increasing cosmological scales and times (redshifts). 
SNIa provide us with a new tool to test how the laws of 
gravity and cosmology were in faraway galaxies [z ~ 0.5). 
The observational limits on G/G come from quite differ- 
ent times and scales [ p8| , |60|j67| ] , but mostly in the local 
and nearby environments at z ~ (solar system, binary 
pulsars, and neutron stars (60| ) . There are also limits de- 
rived from the white dwarf cooling theory ]64| ] , which are 
based on similar arguments to the ones presented in this 
paper. Typical upper bounds give G/G < 10 -11 — 10~ 12 
yr- 1 @. 

Here we have proved by using detailed numerical mod- 
els that if the value of G at the space-time location of dis- 
tant supernovae is different from the local one, it would 
change both the thermonuclear energy release and the 
time scale of the supernova outburst. The change can be 
quantified by means of the change in the Chandrasckhar 
mass M Ch oc G~ 3 / 2 , and our detailed numerical results 
have been interpreted in terms of a very simple physical 
model. To this regard it is important to realize that our 
conclusions would remain unchanged should a modifica- 
tion of the parameters of the explosion lead to a smaller 
mass of 56 Ni synthesized in the supernova event, leading 
to a dimmer supernova. 

We have also shown in a self-consistent way that for 
plausible models for a varying G, such as scalar-tensor 
theories, the possible effect of a varying G on the cosmo- 
logical evolution yields a contribution to the luminosity 
distance relation which is one order of magnitude smaller 
than the effect produced by the same variation of G on 
the Chandrasekhar mass. Thus our approach comple- 
ments the analysis of [ p5p^j6^ | which neglected the ef- 
fects in the physics of the supernovae and used the evi- 
dence for the accelerating universe as a way to constraint 
cosmic evolution in non-standard theories of gravity. 

In this paper we have also found further bounds for 
a varying G from a likelihood analysis of the peak lumi- 
nosities of the Supernova Cosmology Project. Our results 
are summarized in Figures 3 and 4 and Eq.0-|^], with 
values of G/Go < 1.2. We have further translated these 
results into bounds for G/G in Eq.@, G/G in Eq.(l| 
and du>~ 1 /dz in Eq. [jl6|. Some of these bounds are new 
or comparable to other existing estimates from the local 
universe, which typically gives stronger constraints for oj 
or G/Go, at least within JBD models. 

In the context of JBD or STT models the limits we 
find for G/G correspond to oj > 3 — 30 and are there- 
fore less restrictive than the solar system limits oj > 2500 
@. However, STTs could allow for oj = oj{4>). To be 



precise, oj is not required to be a constant, so that oj 
could increase with cosmic time, oj = oj(z), in such a 
way that it could approach the general relativity predic- 
tions {uj — > oo) at present time and still give significant 
deviations at earlier cosmological times [ ^5|^2p3] |. Fur- 
thermore, it has been shown [|55| that the cosmological 
evolution makes STTs practically indistinguishable from 
General Relativity at the present epoch. Our results set 
strong constraints at cosmological distances. 

The interest of these new bounds with respect to the 
other values discussed so far in the literature, is not 
whether or not they are better, but the facts that: i) 
a different method has been tested and used and, ii) our 
bounds correspond to higher redshifts, z ~ 0.5, thus ex- 
tending the constrains on the evolution of G all the way 
from solar/stellar distances to Gpc, that is by more than 
15 orders of magnitude. In this sense, cosmological nu- 
cleosynthesis also offers another limit on the amount of 
variation of G. Generally speaking, the bounds derived 
from primordial nucleosynthesis arise from the sensitivity 
of the abundances of light elements produced at high tem- 
perature to the expansion rate of the Universe at those 
temperatures, especially 4 He. There is a range of opin- 
ions, but there is also the widespread agreement that the 
expansion rate must have been well within a factor of two 
of the standard model. Some might even push for more 
stringent limits that would exclude changes by even as lit- 
tle as ~ 10%, which would be marginally consistent with 
our analysis — see the most recent analysis presented in 
p8[ for a detailed discussion. 

Finally, we would like to stress that new observations of 
distant supernovae, or other standard candles, at higher 
redshifts [z > 1) will constrain even more the current 
limits on the variation of the fundamental constants (see 
Eq. [p^[jl^| ) . To this regard it is important to realize that 
the recently analyzed SNIa 1997ff |p9j , the oldest and 
most distant SNIa ever discovered at z ~ 1.7 [0, could 
provide an important test of the viability of alternative 
theories of gravity. 
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APPENDIX A: SCALAR-TENSOR THEORIES 

The main topic of this paper is how the Hubble dia- 
gram of distant supernovae could help to set constraints 
on a varying G. In order to do that we need to study the 
physics of SNIa, but to be consistent we need to derive a 
luminosity distance relation which also includes the pos- 
sible effects of a varying G on the cosmological evolution 
(sec |7l[] ) . In this Appendix we will show that for a given 
change in G this later effect is typically smaller than the 
one induced by the change in the Chandrasekhar mass in 
the energy release of SNIa. 

The possibility that G could vary in space and/or time 
naturally appears in the framework of Scalar- Tensor the- 
ories of gravity (STTs) such as JBD theory or its ex- 
tensions. These models have recently attracted a large 
interest (see and references therein). To make quan- 
titative predictions we will consider cosmic evolution in 
STTs, where G is derived from a scalar field (f> which is 
characterized by a function uj — uj{<f) that determines 
the strength of the coupling between the scalar field and 
gravity. In the simplest JBD models, uj is just a constant 
and G ~ 4>~ l , however if uj varies then it can increase 
with cosmic time so that uj — uj(z). The Hubble rate H 
in these models is given by: 



8np 
"30" 



1 



a 2 R 2 



3 + 6tf- H 6> (M) 



this equation has to be complemented with the accelera- 
tion equations for a and <f>, and with the equation of state 
for a perfect fluid: p = (7 — l)p and p + 37 Hp = 0. The 
structure of the solutions to this set of equations is quite 
rich and depends crucially on the coupling function u)((f>) 
p7f ). Here we are only interested in the matter domi- 
nated regime: 7=1. In the weak field limit and a flat 
universe the exact solution is given by: 



G = 



4 + 2oj 
3 + 2lj 



r 1 -G (l + z) 1 /( 1+ "). 



(A2) 



In this case we also have that a = (i/t ) (2 " +2)/(3w+4) . 
This solution for the flat universe is recovered in a gen- 
eral case in the limit t — > 00 and also arises as an exact 
solution of Newtonian gravity with a power law G cx t n 
|[72fl . For non-flat models, a(t) is not a simple power-law 
and the solutions get far more complicated. To illustrate 
the effects of a non-flat cosmology we will consi der gen- 
eral solutions that can be parameterized as Eq.|A2| but 



which are not simple power-laws in a(t). In this case, it is 
easy to check that the new Hubble law given by Eq.[Al] 
becomes: 



H 2 = H 2 



Cl M (l + z) 3+1 /( 1+ ") + f2 fl (l + zf + Q A 



(A3) 

where £Im,Qr and JIa follow the usual relation: J7m + 
fin + (Ia = 1 (an overall factor would just redefine the 
value of Ho) and are related to the familiar local ratios 
(z -» 0): f> M = 87rG p /(3ff 2 ), Sl R = 1/{RH ) 2 and 
n A = A/(3H 2 ) by: 



Om = 



M 
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3 + 2uj 

4 + 2w 

1 1 



(l+Uj) 6 (l+Uj) 2 



;A A =2i;ft JI =fk (A4) 
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(A5) 



Thus the general relativity limit is recovered as uj — > 
00. For a Hat universe, the luminosity distance dr, — 
cIl(z,CImiQa,w) is related the (line-of-sight) comoving 
coordinate distance r(z) as: 



r{z) 



H 



■(1 



dz' 



(A6) 



In the general case we have to replace the integral with 
its trigonometric or the hyperbolic sinus to account for 
curvature |7q| . In the limit of small z we recover the usual 
Hubble relation: d Hq/c = z — (1 + qo)z 2 /2 where a new 
deceleration go parameter is related to the standard one 
by: 



qo , 



(A7) 



One can see from these equations that even for relative 
small values of uj the effect of a varying G on 4 is small. 
For example for the flat case (Or = and CIm — 1) at 
z = 0.5 we have c?l = 0.5456 for uj — 10 and dz, = 0.5505 
for uj — 10 3 (g?l in units of c/Hq). Thus, the change 
in uj produces brighter apparent objects, in this case 
Am = —0.019, which would tend to partially compensate 
the dimmering produce by the a varying G on the Chan- 
drasekhar mass: in this case Am — +0.060. In general, 
we find that the cosmological effect in the Hubble dia- 
gram of SNIa is always smaller, by factors of a few, than 
the effect produced by a varying G on the Chandrasekhar 
mass. Also in the general case, the cosmological evolution 
in a model with increasing G at high z tends to decrease 
the acceleration (with respect to the case with constant 
G), which partially compensates the apparent increase 
due to the the change in the Chandrasekhar mass. 
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